On the thermodynamics of first-order phase transition smeared by frozen disorder. 
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The simplified model of first-order transition in a media with frozen long-range transition- 
temperature disorder is considered. It exhibits the smearing of the transition due to appearance 
of the intermediate inhomogeneous phase with thermodynamics described by the ground state of 
the short-range random-field Ising model. Thus the model correctly reproduce the persistence of 
first-order transition only in dimensions d > 2, which is found in more realistic models. It also allows 
to estimate the behavior of thermodynamic parameters near the boundaries of the inhomogeneous 
phase. 
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The smearing of first-order phase transitions by frozen 
disorder of random bonds (random transition tempera- 
ture) type is ubiquitous phenomenon, which can be ob- 
served in crystalline solid solutions, porous media, gels 
and composites. The jumps of thermodynamic parame- 
ters proper to such transitions could be diminished or 
completely wiped out by the disorder. According to 
heuristic criterion derived in Ref. Q this takes place for 
sufficiently large range of disorder correlations, greater 
than the order parameter correlation length and the 
length defined by the ratio of interphase surface tension 
and latent heat. The rigorous result was obtained in Ref. 
0for random bond Potts models (RBPM) with first-order 
transitions in space dimension d = 2. It was established 
that latent heat vanishes in planar Potts models irrespec- 
tive of disorder strength while long-range order persists 
at low temperatures. Further numerical studies of q-state 
]^gp]y[3j4j5j6 ];ja,ve shown that in d = 3 first-order transi- 
tions are wiped out only at sufficiently large disorder. 

This situation can be compared with that in random 
field Ising model (RFIM), where, according to Imry-Ma 
arguments^ and rigorous resultsSiS, transition can also ex- 
ist at c? > 2 only. This hints on the possibility of RBPM 
first-order transitions and RFIM to belong, in some sense, 
to the same universality class. This possibility is strongly 
corroborated by the equivalence of RBPM at g > > 1 and 
RFIM established in Ref. |^ for d — 2. Qualitative argu- 
ments in favor of close relation between first-order transi- 
tions in random bond systems and RFIM were advanced 
in Ref. 0. Yet to date there are no analytical results 
on the thermodynamics of generic smeared transitions, 
revealing such relation. So one may try to obtain some 
insight in this problem considering simplified models. 

Here we study the thermodynamics of simple model 
with long-range correlated disorder, which may capture 
the main features of the phenomenologjsi. According to 
Ref. T the smearing in random bond (random transition 
temperature) system results from the appearance of in- 
homogeneous equilibrium state, consisting of clusters of 
two phases (ordered and disordered), which is energeti- 
cally more favorable than the homogeneous one. We find 
that in present model there is a definite temperature in- 
terval where such intermediate inhomogeneous phase ex- 



ists and its thermodynamics is described by the RFIM 
ground state for the strong first-order transition. We ob- 
tain the estimates for thermodynamic parameters near 
the boundaries of this phase. 

Let us consider a d-dimensional sample undergoing 
first-order phase transition and divide it on hypercubes 
of size L in lattice units. We assume that in each hyper- 
cube the transition temperature is with probability p 
and T_ with probability 1 — p, T_ < r_|_. If L is much 
larger than the order parameter correlation lengths in the 
ordered and disordered phases, 



L>max(a,6) > 1, 



(1) 



thermodynamics of each cube can be described by 
the density of inequilibrium thermodynamic potential 
f{ip,T±) of infinite sample. Here y is multicomponent 
order parameter and t± = T/T± — 1 denotes the reduced 
temperatures, which correspond to the regions with tran- 
sitions at r_|_ and T_ . So the effective Hamiltonian of 
the system can be expressed as 



JV 



H iW} , a) = L'^ ^ / (^„ T,. ) + Wi„t iW}) . 



(2) 



with random variables ai = ±, N being a number of 
cubic blocks and Hint {{f}) representing the surface in- 
teraction of neighboring blocks. The interaction term is 
proportional to L'^^^ and it should depend on the differ- 
ence of if in the nearby cubes. So we put 



Hint {{v}) = L 



d- 



(3) 



where A is some constant. 

The density of the average equilibrium potential of the 
model is 



N-^L-''T (in / dtp exp[-/3H ({¥>}, cr) 



(/) 

(4) 

As L ^ 1 integral in Eq. Q can be estimated via the 
method of steepest descend. So we must find a global 
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minimum of 7i ({y'} , c). Here we assume that at temper- 
atures including the interval T_ < T < both / {ip, r_|_) 
and / (c^, r_) have minimum at if — and q minima at 
<fi ^ <Pa (t+) or (f = ifa (t_), a = 1, 2, q, correspond- 
ingly. We also suppose, that / (ip, r) is invariant under 
some symmetry group, so for all a 

fa (t) = LPs (r) , f{ip^,T)^ fs (r) 

and J2 fa = 0. 

Then H {{ip} , cr) has {q + 1)" local minima at which 
every ipi is close to one of the f {ip,T„^) minima differing 
from them by terms of order 1 / L. These local minima can 
be parameterized by the vectors i/^i, which take values 
ifa/'^s and 0, so we have at them 

ip,^ip,ips (r^J-f 0(l/i). 

For the strong first order transition we may assume 
that at T_ < r < T+ 

ips (t+) W ips (t_) = ips. 

Then we get for the value of Hamiltonian at local mini- 
mum 

N 

1=1 

+2i''-M^2^(V^,-t/,,t/>,) (5) 

Here /o (r^) = / (0, r^), Sf (r<^) = fs (r^) - /o (r^) and 
we used the relation ipf = ipi being the consequence of 
V^i = {0, 1}. So we have 

{f)=N-'L-''(mmH{m,a)^ 

Apparently, the global minimum is realized on the con- 
figurations, which have parallel ipi and tpj on neigh- 
boring sites (if they are both nonzero). Limiting the 
choice to such configurations, we can put tpiXpj = ipi'4'j 
in Eq. |(SJl. Then, expressing ipi via Ising spins Si = ±1, 
ijji = (1 - Si) /2, we get 

{f)^l{{foM), + {fsM)J+dJ 

+N'^ (^mmE{s,(T)'^ , (6) 
^ J 

E{s,(t) = -^i?iS, - ^^SjSj, 

H,^Sf,j2, J = A^l/L. 

Thus for large disorder correlation length L and strong 
first-order transition the average equilibrium potential of 



this simple random-temperature system is determined by 
the ground state of short-range Ising model with random 
fields having bimodal distribution. 

The present results rely on the large disorder corre- 
lation length L and the possibility to ignore the order 
parameter temperature dependence for strong first-order 
transitions. Actually it could be sufficient to require 
(along with the condition the first-order transition to 
be strong with small order parameter fiuctuations around 
the local minima. Indeed, Eq. © holds, when 

L-^ EE min (0, 0) , {VsM » (7) 

where /^j„ t) is the lowest eigenvalue of the matrix 
of the second-order c^-derivatives of the potential. So, if 
Lc < max (^1,^2) (meaning small order parameter fiuc- 
tuations), the condition Q does not impose further re- 
striction on L as compared with that of Eq. Also 
small Lc implies only slight temperature dependence of 

{t±)- In case of a soft transition, Lc S> max(^i,^2), 
and for max (^1,^2) ^ L < Lc one should take into ac- 
count the fiuctuations around the local minima, as well 
as ips (''"i) temperature dependence, which would essen- 
tially modify the effective spin Hamiltonian breaking the 
simple RFIM equivalence. 

To find the RFIM ground state in Eq. © IS non- 
trivial task only in the interval T_ < T < owing to 
the definite temperature dependence of two possible Hi 
values. Indeed, H+ = Sf+ll and = (S/-/2 grow 
monotonously with increasing T and change sign at T4. 
and r_ correspondingly. Particularly at small \t± \ « 1 
we have 

6f± ~ T±T±dS, 

where SS > is the entropy jump at the transition in 
either component. So at T > all Hi > and ground 
state is Si — +1 for all i and 

(/> = (/o(r.)), = p/o(r+) + (1 - p) Mr-), 

while at T < r_ all Hi < , so ground state is Si = —1 
and 

(/) = (./s(t.)),. 

Thus above T+ all cubic blocks are in the disordered 
phase and equilibrium potential is that of this phase av- 
eraged over random transition temperature. Similarly 
below T_ we have all blocks in the ordered state. So in 
the interval T_ < T < T_|_ gradual transformation from 
homogeneous disordered state to a homogeneous ordered 
one takes place. Inside this interval the remnant jumps 
of first-order transition can exist at some Tq defined by 
the condition (Hi)^ = or 

(/o(T-.)), = (/.(r.)),. 

The present results can be trivially generalized to the 
case when random transition temperatures are continu- 
ously distributed between T_ and r+, which would give 
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continuous Hi distribution with bounded support. Yet, 
according to numerical studiesiitiSiiiii, thermodynam- 
ics of RFIM at T = does not crucially depend on the 
random field distribution, apart from the macroscopic 
ground state degeneracy in the bimodal case-^. As gen- 
erally the long-range order in the RFIM ground state ex- 
ists only in d > 3 (for small disorder), we have the same 
condition for the persistence of first-order transition in 
the present model for generic transition temperature dis- 
order. So it is similar in this respect to the more realistic 

j^g p ]y[2j3j4j5j6j9jl0jl6 

Yet there is apparent qualitative discrepancy with nu- 
merical results for RBPM, which indicate that \n d — 
2 first-order transition transforms into a second-order 
one2*i2ii&, while the study of 2d RFIM ground state in 
homogeneous field^^ shows that there is no transition at 
{Hi)^ = [T — Tq). Thus present model cannot eluci- 
date the origin of the instability and the nature of the 
order parameter appearing at To in 2d RBPM. Still we 
may suppose that far from Tq RFIM ground state cor- 
rectly reproduce the qualitative features of the inhomoge- 
neous equilibrium state in realistic models, especially the 
existence of intermediate inhomogeneous phase in defi- 
nite temperature interval T_ < T < T_|_ in which random 
transition temperatures vary. Actually, the Imry-Wortis 
phenomenologyi already implies its existence, but to date 
numerical2i^*^*S*i& and renormalization group2ii2iii stud- 
ies of RBPM were not intended to reveal such phase. 

We can give rough estimates for the behavior of ther- 
modynamic parameters near T_ and inside the inho- 
mogeneous phase. Thus near r+ negative fields = 
(5/-f/2 are small so in the sea of positive spins only very 
large (and, hence, very rare) clusters of negative spins 
(ordered phase) can have sufficiently large energy gain as 
compared to the surface energy loss to diminish RFIM 
ground state energy. Then we can choose approximate 
trial ground state dividing the lattice on the cubic blocks 
of size A >> 1 and putting all spins positive except for 
the cubes in which homogeneous component of random 
field He = A^'^ Y.Hi < ~dJ/A. In the last blocks we 

put all Si = —1. As there is small probability to find 
nearby two or more such blocks, we get for the energy 
of such spin configurations in the most of the disorder 
realizations 



N 



E (A, a) = -NdJ - ^ H, 



N/A" 



-A'' J2 + dJ/A) d {-He - dJ/A). (8) 



As 



(mmE{s,(T)) <(mmE(A,a)) < mm (E (A, a)) ^ . 

\ s /(t\A / a a 



culating nnn (A, cr))^. From Eq. (jSJ it follows 



mm{EiA,a)}jN 



U+ = min 

A 



-dJ - (H,)^ + U+ 
J dHW (H) {2H + dJ/A) 



(9) 

(10) 



Here W [H] is the distribution function for homogeneous 
field He = A^"* ^ Hi in blocks considered. For H close 

to H+ and A >> 1 it has the form 



W{H) = 



A" 



2tt {H_-H+){H-H+) 



1/2 



P 



The minimum of the integral in Eq. H10|l is attained at 
the size 



A, 



dJ 

W7\ 



1 



\H^ 



2d\n{l/p) {djy 



which diverges when T — > T-i_ — 0. Thus we have near 

\ 1/2 



7rdln(l/p)fl"_ 



1/2 



47rdln(l/p) 5f- 



exp 



ln(l/p) 



2dJ 



(11) 



This expression describes, apparently, the contribution 
of very large and very rare clusters of ordered phase ap- 
pearing immediately below r_|_, which diminish the en- 
ergy. The probability to find such clusters, , vanishes 
very fast at T_|_ (as exp ^— const/ \tj^ j'^^ ) and it dominates 

the temperature behavior of U+. On physical grounds, 
the presence of this dominating term in Eq. Hll|) can be 
expected also in the true ground state energy, so we may 
assume that expression in Eq. © is rather close to it. 
So we suppose that near r_|_ 

(/)«(/i(r.)), + ;7+. 

Then average entropy and heat capacity are 

d+l 



{S) « {So (r.)). 



2 J \\Sl 



2dJ 



{C) « (Co (r.)), + 
Here So (r) = - 



2dJ 



\U+\T^ 

4.P y\5i 

dfo (r) 
dT ' 



2d+2 



Coir) 



5S\n{l/p). 

5S^ \n{l/p) 
.d'fo (r) 



-T- 



QT2 



Thus slight diminishing of the entropy starts already at 
T_|_ indicating the rounding of the transition. Accord- 
ingly, the heat capacity becomes larger than its mean- 



we can get upper bound for the ground state energy cal- field value below T+. 



4 



In order to find the spontaneous order parameter, one 
should consider a system in an infinitesimal field h con- 
jugate with the order parameter. In the presence of such 
field wc should take into account only one nontrivial min- 
imum of the potential, namely, the one with ipa direction 
closest to that of h as the others have higher values. Let 
it be the minimum at ipi. Then two values of spins in Eq. 
© correspond to <^ = (s = -1-1) and (p = (fi {s = —1), 
so the average spontaneous order parameter is 



Sf3 ^ 1/2 



<P1 



Here s^^"* (cr) is ground state spin configuration. When 
there are several ground states the average over all of 
them should be taken in this expression. Using the above 
trial ground state with A = A+ we get 



(if) = 



dHW [H] « ipi \U+/H+\ 



(12) 



Thus spontaneous order parameter appears continuously 
below T-|_. Yet for every disorder configuration the transi- 
tion into inhomogeneous phase near this point is a first- 
order one. The absence of discontinuities at T+ in the 
average thermodynamic parameters results from the van- 
ishing probability to have the finite jumps in all of them 
including {ip). 

In the vicinity of T_ , where rare clusters of disordered 
phase exist, analogous treatment gives the following es- 
timates 



U. 



{if) « <pi (1 - C/-/i/-) 



d-l-1 



(55 In 



1-p 



47rdln(l -p) 5f^ 



exp 



Thus near T_ there also are exponentially small contri- 
butions (proportional to exp {— const /t'^^) to the mean- 
field values of thermodynamic parameters, which indicate 
the softening of the transition at Tq. 

We should note that the slight anomalies at the bound- 
aries of inhomogeneous phase are hard to observe in the 
numerical simulations. As they result from the rare ap- 
pearance of very large clusters of opposite phase, the 
large samples and large number of disorder realizations 
are needed to reveal these anomalies at T+ and T_. Yet 
the tails in the temperature dependencies of thermody- 
namic parameters, similar to those described by Eq. H12|l . 
are often seen in experimental studies of first-order tran- 
sitions in random media (see, for example, Ref. ilStl . This 
may indicate the presence of intermediate inhomogeneous 
phase in real systems. 

This may also imply that the present model with seem- 
ingly unrealistic disorder is more closely related to the 
realistic models than one may expect. Indeed, one can 
imagine that application to, say, RBPM of some sort of 
renormalization group procedure, which eliminates the 
order parameter fluctuations on scales smaller than some 
large L, will result in effective Hamiltonian, similar to 
that in Eq. ^ in some range of the model parameters. 
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